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, Abstract 



We study the asymptotic behaviour of the Bianchi type VIo universes 
with a tilted 7-law perfect fluid. The late-time attractors are found for 
' the full 7-dimensional state space and for several interesting invariant 

subspaces. In particular, it is found that for the particular value of the 
equation of state parameter, 7 = 6/5, there exists a bifurcation line which 
signals a transition of stability between a non-tilted equilibrium point to 
5-H . an extremely tilted equilibrium point. The initial singular regime is also 

discussed and we argue that the initial behaviour is chaotic for 7 < 2. 



X ■ 1 Introduction 



In a recent paper we analysed how tilted fluids affect the general late-time be- 
haviour of cosmological models of Bianchi type [1] . As the asymptotic behaviour 
of all Bianchi models with a non-tilted 7-law perfect fluid are determined [2-4] , 
the sensitivity of the non-tilted equilibrium points was checked with regards to 
inclusion of tilted fluids. In this paper, we will consider one of the models in full 
generality; namely the Bianchi type VIq model containing a tilted 7-law perfect 
fluid. 

The last decades, there have been some investigations of Bianchi models 
with a tilted fluid [5-7]; in particular, type II [8] and type V [9-11]. Apart 
from these special Bianchi types, the general behaviour of tilted Bianchi uni- 
verses seems to be poorly understood. The work [1] gave us some hints where 
interesting behaviour may occur, but a more elaborate analysis is needed in 
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order to understand the more general behaviour of universes with tilted flu- 
ids. More specifically, in [1] it was shown that for the Bianchi type VIq model 
there are no stable non-tilted equilibrium points for equation of state parameter 
obeying 7 > 6/5. The value 7 = 6/5 signals the onset of an instability with 
respect to tilt implying that at late times the peculiar velocities of the fluid is 
a major contributor to the cosmological shear. The Bianchi type VIq is not the 
most general model but is sufficiently general to account for many interesting 
phenomena. Thus we will assume that the universe contains a 7-law perfect 
fluid (0 < 7 < 2) which in general can be tilted. The state space is thus 7- 
dimensional [2] (compared to the 8-dimensional state space for the most general 
models) . 

The dynamical systems approach - which is the method adopted here - 
requires that we find all equilibrium points for the system of equations. These 
equilibrium points play a special role in the evolution of the system. Not only 
do they correspond to exact self-similar solutions to Einstein's equations, but 
some of them may even serve as attractors for more general solutions. However, 
not many exact solutions with tilted fluids are known for the Bianchi models. 
Some type II solutions have been found [12], and Rosquist and Jantzen found 
some type VIq solutions [13, 14]. Here, we will show that the solutions (for 6/5 < 
7 < 3/2) found by Rosquist and Jantzen are attractors in a four-dimensional 
invariant subspace of the full state space. However, in the full state space, there 
always exists an unstable mode. In the full state space, the future attractor for 
6/5 < 7 is an extremely tilted model which is connected via a line bifurcation 
at 7 = 6/5 to a non-tilted equilibrium point. This 7 = 6/5 bifurcation, which 
correspond to a one-parameter family of exact tilted solutions, seems to have 
been first discussed in a recent paper by Apostolopoulos [15]. However, not 
much discussion is devoted to this solution as the solution itself is implicitly 
given in terms of a cubic. Here, we have introduced a different parameterisation 
which makes it possible to explicitly write down the solution in terms of the 
expansion-normalised variables. 

We also discuss the initial singular regime, which is a lot more complicated. 
In fact, there are indications that the initial singular regime possesses an oscilla- 
tory, very likely even a chaotic behaviour similar to that found in other Bianchi 
models [16-24]. 

The paper is organised as follows. Next, in section |21 we write down the 
equations of motion using the orthonormal frame formalism. In sectional we 
discuss the equilibrium points which are important for the late-time asymptotic 
behaviour. Then, in sectional we discuss their stability and present some general 
results regarding the asymptotic behaviour of the tilted type VIq model. The 
initial singular regime is discussed in section [S] before we finally summarise our 
results in sectional 
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2 Equations of motion 

The energy-momentum tensor for a tilted perfect fluid is 

T^v = {P + 'p)UnUu + PQiiv, (1) 

where u'^ = (cosh/3, sinh/3c") is the fluid velocity. The spatial vector c" is 
chosen to be a unit vector in the tangent space of the surfaces of homogeneity; 
i.e. c"ca = 1. We will further assume that the fluid obey the barotropic equation 
of state, 

p = (7 - < 7 < 2. (2) 

In terms of the unit normal vector u = eg to the group orbits the energy- 
momentum tensor takes the imperfect fluid form 



where 



T^tu = (p + p)u^u^ + pg^y + 2g'(^u^) + tt^^, (3) 



p = (l + 7sinh2/3)p, (4) 
p = ^7-1 + ^7 sinh^ p, (5) 
Qa = 7PCOsh/3sinh/3ca, (6) 



TTab = J p sinh /3 yCaCb- -habj ■ (7) 

For the Bianchi cosmologies - which admit a simply transitive symmetry 
group acting on the spatial hypersurfaces - we can always write the line-element 
as 

where a;" is a triad of one-forms obeying 

1 

-( 
2 

and C°j^ depend only on time and are the structure constants of the Bianchi 
group type under consideration. The structure constants can be split into 
a vector part at, and a trace- free part n"*" by [25] 

= Sbcdn'^'' - 5%ac + S^a^. 

The matrix rf'^ is symmetric, and, using the Jacobi identity, = {1/2)C\^ is 
in the kernel of n"^ 

n"''a6 = 0. 

For the type VIq model, Oc = and Uab has two non-zero eigenvalues with 
opposite sign. This implies that we can choose a frame such that the structure 
constants can be written 



nab 




n n 
n fi 



ab = 0. (8) 
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Furthermore, the type VIq has < n^. ^ The equations of motion can now be 
written down. 

Fohowing [2] we introduce expansion-normaUsed variables to write the sys- 
tem as an autonomous system of differential equations. In the notation of 
ref. [26] we define 



— 2E+ \/3Si2 



\/3I]i3 

'^ab = ■\/3Ei2 S+ + a/SE- V^^23 

_\/3Si3 \/3S23 — \/3E_ 

N22 = -/V33 = VSN, N23 = N32 = VSN. 



We also introduce the three-velocity V by 

V 

sinh/3 = 



< y < 1. 



For the expansion-normalised variables the equations of motion are: 
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T' — 



-•IS 
^23 

N' 
N' 

n' 
v 



(<?-2)S+ + 3(S?2 + I]?3)-27V2 

-t-i [ANT,-vi + {NY.12 + iVSis)^ - (A^Sis + ^"212)1^3] 
(g - 2)S_ + V3(E22 - E^a) - 2iiiS23 

+^ [(ArSi2 + A^Si3)t;2 + (A^Si3 + Ni:^2)v3\ 
3E+ - \/3E_ -|- \/3iVwi) E12 

- + V3E23 - V3iVi;i) Sis 
(q-2- 3S+ -h \/3E_ - y/SNvi^ S13 

- + V3Y.23 + VSNvi'^ S12 

{q - 2)S23 - 2\/37ViV + 2i?iS_ 2\/3Si2Si3 
-\/3u2(Ari]i3 + NT, 12) 
{q + 2S+) N + 2V3S23^ 
(9 + 2S+) TV -h 2v/3S23Ar 



(9) 



(10) 



1 + (7 - 
V{1 - v^) 

1 - (7 - 



(11) 

(12) 
(13) 

(14) 

(15) 
(16) 
(17) 

[2q - (37 - 2) + [2(z(7 - 1) - (2 - 7) - 7^] F^} (18) 

2 [(37 - 4) - 5] (19) 



Note that the Bianchi type II is the limit where = n^, and that Bianchi type VIIo has 
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where 



S2 




s 



Ea b a 1 a T r a 

abC C , C Ca = I, V = VC 

vl+vi+vl 



(20) 



These variables are subject to the constraints 
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2I]_iV[l + (7-l)y2] +jQvi 

- (SiaTV + Si3^) [l + (7 - 1)1^'] + 1^V2 

(EigA^ + Sl2iV) [1 + (7 - 1)1^'] + l^Vs 

NRi - VS^^N. 



(21) 
(22) 
(23) 
(24) 
(25) 



The parameter 7 wiU be assumed to be in the interval 7 e (0, 2). 

Note that Ri, which is the component of the rotation tensor describing 
rotations with respect to the axis ei, is implicitly defined via ea. H25|) . For 
the type VIq model we will in practice solve this equation by introducing the 
parameter A instead of TV, by TV = XN. For the type VIq model, this parameter 
is bounded by < 1. 

2.1 The state space 

The constraint (|21|l implies that S±, S12, S13, S23, N and H. are all bounded. 
Combining the equations for TV and TV we get the equations 



Thus if the initial data have TV^ < TV^, then this will hold for all times. 
Hence, for type VIq, TV will be bounded as well. The invariant subspaces 
TV ± TV = correspond to Bianchi type II universes. Note that, given the 
set (Ei, E12, Si3, S23, TV, TV) we can determine from eq. and v^, 

and from eqs. (|22|l . H23(l and (|24(l . respectively. We also require < < 1 
for physical reasons (we are not allowing superluminal velocities). This implies 
the bounds 



However, as can be shown, the last of these inequalities is redundant due to the 
fact that the first two imply the third (see Appendix Hence, the state space 
can be considered a subspace of a compact region in K.^. 




(26) 



+ + + Sjg + S^g + TV^ < 1 

TV2 < TV^ 

4E2^TV2 + (l]i2TV + Ei3TV)% (Si3TV + I]i2TV)^ < 1. 



(27) 
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There arc also some discrete symmetries under which this system is invariant. 
These are (of course, also compositions of these maps are symmetries) 



(pi : (S_, Ei2, Sis) 
02 : (Si2, S23, N) 1 

4>3 ■ {^13,^23, N) 1 




(28) 
(29) 
(30) 



To understand the physical importance of these maps we can apply them to the 
vector u": 



Hence, these discrete symmetries correspond to different permutation and in- 
versions of the axes. These symmetries are intrinsic to the type VIq geometry. 

2.2 Invariant subspaces 

The invariant subspaces play an important role in the evolution of Bianchi type 
VIq universes. They are important to study since they are invariant under the 
evolution of the system; i.e. if a point p lie in one invariant subspace, then also 
the maximal extended evolution of p will lie entirely inside it. 

Some of the physically interesting invariant subspaces are as follows (we will 
also assume that the boundaries are included): 

1. T{VIoy. The full state space of tilted type VIq. 

2. F{VIo): The set of fixed points of the map <pi. Given by S_ = 0, S12 = 



3. T2{VIo): A Bianchi type VIq with a two-component tilted fluid. It is 
the set of fixed points of (^2 (or (j)^). Given by S12 = S23 = N = {or 
Sis = ^23= N = 0). 

4. Ti(y/o): Bianchi type VIq with a one-component tilted fluid. Given by 
S12 = Sis = 0. This is the fixed-point-set of 02 o (/'s- 

5. B{VIo): Non-tilted Bianchi type VIq. Given by S_ = S12 = Sis = V = 
0. 

6. T^{II): The tilted type II boundary. Given by iV = ±iV. Note that these 
two subspaces can be mapped onto each other using 02 or 03. Because of 
this, we will in most cases not differentiate between these two invariant 

subspaces. 

7. B{II): Non-tilted type II. Given by = and V = 0. 

8. B{I): Bianchi type I universes. Given hy N = N = V = 0. 



01 

02 
03 



{V1,V2,V3) 
{V1,V2,V3) 
{V1,V2,V3) 



{-Vl,V3,V2) 
{Vl, -V2,V3) 
{vi,V2,-V3) . 



-E 



13- 
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Subspace 


Dim 


Subspace 


Dim 


TiVIo) 


7 






F{Vh) 


5 




6* 


T2{Vh) 


4 


Bill) 


4* 


TiiVIo) 


5 


B{I) 


5* 


B{VIo) 


4 


dvB{I) 


4* 



Table 1: The dimension of the state space for different invariant subspaces. The 
asterisk indicates that the true number of physical degrees of freedom should 
be reduced due to unphysical gauge freedoms. 

9. dvB{I): The Bianchi type I vacuum boundary. Given hy N = N = V = 

r2 = 0. 

3 Equilibrium points of importance for the late- 
time behaviour 

The system of equations has numerous equilibrium points which correspond to 
fixed points in the evolution of the system. That is, if the state variables are 
written as a vector X, and we write the system of equations as 

X' = F(X), (31) 

then a point Xq is an equilibrium point if 

F(Xo) = 0. 

These equilibrium points usually play a particular role in the evolution of the 
system. Not only are they exact solutions, but some of these equilibrium points 
may act as past or future attractors for the set of solutions. Hence, for example, 
at late times these equilibrium points may be a good approximation of more 
general solutions. 

Some of the physically interesting equilibrium points for the future evolution 
are given below. If some equilibrium points are related via one of the discrete 
symmetries, permutations, or other gauge symmetries, we only give one of them 
because the solutions they represent are equivalent. We also indicate the small- 
est invariant subspace mentioned in section [2. 21 to which they belong. 

3.1 Non-tilted 

1. X{I): FRW_ 

^ N ^ N i(37 - 2). 

Invariant subspace: B{I). 
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2. CS{II): Collins-Stewart type II (2/3 < 7 < 2) 

S_ = S12 = Si3 = y = 0, E+ = - j^(37 - 2), iV = AT 

S23 = -^(37 - 2), = ^(37 - 2)(2 - 7), 

n=^(6-7),g=i(37-2). 

Invariant subspace: B{II) 

3. C(F7o): Collins VIq (2/3 <_ 7 < 2) 

I]_ = E12 = S13 = E23 = iV = y = 0, S+ = -i(37 - 2), 

= ^(37 - 2)(2 -j),n= 1(2 - 7), « = 1(37 - 2) 
Invariant subspace: B{VIo) 

3.2 Intermediately tilted 

1. n{II): Hewitt's tilted type II (10/7 < 7 < 2) [12] 
S_ = 0, S+ = 1(97 - 14), S23 = -^(5 7 - 6), 

y,„_y___V6 / (2-7)(ll7-10)(77-10) /u - 

^12 - ^13 - — §-V vFTTs ' -'^ - -'Vj 

Ar2 _ 3(2-7)(57-4)(37-4) , .2 _ (37-4) (77- 10) 
~ 4(177-18) ' ^ ~ (ll7-10)(57-4) ' 

" = 4(!^(2l7' - 247 + 4), g = i(37 - 2). 
Invariant subspace: T^{II) 

2. £(//): Type II line bifurcation (7 = 14/9) [8] 

S± = 0, E23 = - ^,N^N, ^ j^{2b + 1)(17 - 86) 

Ei3 + S12 = -1^^(46 +1)(8- 36), Si3 - E12 = -f^6, 

V' = g7-t6i!s-V.! ^ ^ = ill (16^^ - + 59), 
0<6<1, g=|. 
Invariant subspace: T^{II) 

3. 7^(^//o): Rosquist-Jantzen (6/5 < 7 < 3/2) [14] 

Ei3 = E23 = ^ = 0, E+ = -i(37 - 2), g = i(37 - 2), 

1 [ (57 - 6) [97^ - 137 + 6 + 3(2 - 7)5] \ 
~ 2| 7(97-10) j 

1 f (2 - 7) [4572 - 657 + 18 - 3(57 - 6)s] 
~ 2 I 7(97 - 10) 

C3 = 

where s = \/{') — 1)(97 — 1). 

E_ = —-5c\, E12 = 5c2 
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where S = \/3[(37 - 2) - 2(2 - j)4]/{6c-i4). 

y = A/TV, = 1 - E - iV^, 
where N <0, and 



Invariant subspace: T2(F/o) 

4. /:(F/o): Type VIq line bifurcation (7 = 6/5)^ 

-|, ci=0, g=|, 



S_ - E23 - 0, S+ = -f , ci = 0, g = 4 



V _ _v — J_ / l8V^^-6A(5V'^+4) Q _ 3 (V^+5)(1+5A) 
^12— ^13— 10 V 2T73+3A+i ' " ~ 25 2y^+3A+l ' 

^2 _ 15 + 211A + 397A2 + 25A3 ± 3(5A + l)^/s 



10(3-5A)(1 - A)2 



where 



(A + 24 - 9V5f - 16(61 - 27\/5)] 
X (A + 24 + 9^5)2 - 16(61 + 27^5)] 

-24-9\/5 + 4-\/61 + 27\/5 < A < 0. 
Invariant subspace: F{VIo) 

3.3 Extremely tilted 

1. £(//): Extremely tilted type II (0 < 7 < 2) 

V , — n V — 2\/3 V _1_ V — 10 
2j± — U, 2^23 — q-j ^13 + ^12 — ^^\/ 



Si3 - E12 = -l^f, = TV = y = 1, = ^ ^ 4_ 
Invariant subspace: T+(//) 



2. £1(^/0): Extremely tilted type VIq (0 < 7 < 2) 

s/I V _i AT — 

2' 



y = 1, Si3 = S23 = TV = 0, s+ = -|, s_ = ^, S12 = -i, TV = -i 

Invariant subspace: 72(^/0) 

^Thcsc solutions seem to be found first by Apostolopoulos [15]. However, he writes the 
solutions implicitly as solutions to a cubic. As can seen, we managed to write them explicitly 
down in terms of one free variable which makes it easier to interpret the solutions. 
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Figure 1: The web of future stable equilibrium points. 




Invariant subspace: F{VIo) 



4 The late-time asymptotic behaviour 

An important key to understanding the late-time behaviour of the dynamical 
system is to consider the future stable equilibrium points. Assume that Xq is an 
equilibrium point; i.e. F(Xo) = 0. Then we can write the equations of motion 
as 

SX.' ^ ASX + OiSX.'^). (32) 

The local stability of the equilibrium point, Xq, depends on the eigenvalues 
of A; if all eigenvalues have negative real parts, then Xq is a future attractor. 
Every eigenvalue with a positive real part signals an unstable mode. 

The (local) future stable equilibrium points are for the different invariant 
subspaces given in Table |3 and the various connections between them are il- 
lustrated in FigH The eigenvalues of the linearised system are discussed in 
Appendix IbI 

For the non-tilted subspaces B{I), B{II), and B{VIq), it is proven that 
generic solutions approach the respective local future attractors [2]. In the 
tilted case we can only partially prove that the local attractors are also global 
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Space 


Matter 


Stable point 






T(T\ 


ft(TT\ 




T(T) 








1 (11) 




X[l) 










10/7 < 7 < 14/9 


nil) 




7 = 14/9 


£{11) 




14/9 < 7 < 2 


£{II) 


B{VIo) 


< 7 < 2/3 


I{I) 




2/3 < 7 < 2 


CiVIo) 


TiiVIo) 


< 7 < 2/3 


1(1) 




2/3 < - < 2 


C(VIo) 


T2{VIo) 


< 7 < 2/3 


1(1) 




2/3 < 7 < 6/5 


C{VIo) 




6/5 < 7 < 3/2 


n{vio) 




3/2 < 7 < 2 


SliVIo) 


F{VIo) 


< 7 < 2/3 


I{I) 




2/3 < 7 < 6/5 


C{VIo) 




7 = 6/5 


CiVIo) 




6/5 < 7 < 2 


£2{VIo) 


T{VIo) 


< 7 < 2/3 


I{I) 




2/3 < 7 < 6/5 


C{VIo) 




7 = 6/5 


C{VIo) 




6/5 < 7 < 2 


SiiVIo) 



Table 2: The future stable equilibrium points for various invariant (sub)spaces. 
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attractors. The subspace T^{II) was discussed in [8]. We will use these results 
in the following way. Assume we have a subspace, S (given by y = 0, say), for 
which a function Z is strictly monotonically decreasing; i.e. 

Away from 5, we can write Z' = a{y)Z where a(0) < e < 0. The function a{y) 
is continuous in y; hence, there will exist a, S > such that 

a(y)<0, for \y\ < S. (33) 

So suppose the solution curves have the property y 0, then after sufhcently 
long time the function Z will be monotonic along the solution curves. The late 
time analysis can therefore to some extent be extracted from the monotonic 
functions in S. 

In the following, uj{p) will mean the future asymptotic set; i.e. for a solution 
curve c(t) having c(0) — p, then Lu{p) = lim^^oo c(t). 

Theorem 1 (No-hair). For 7 obeying < 7 < 2/3, any p e T{VIo) for which 
V <1, andn>0 has uo{p) = I{I). 

Proof. Part of the proof relies on the observation that the function S = E^bC^c^ 
obeys the bound \S\ < 2E < 2 (see Appendix 0. From the equation for V we 
thus have 

Hence, if < 7 < 2/3, V will be monotonically decreasing. 

There is also an other monotonically increasing function, namely 

{(3ny = [2q-{3^-2)]{f3n), (34) 

(1 - ^2^(2-7) 



Note that 



^2 , ,0 o.A^r2 , 7(4-37)1^^ 

i + (7-i)y2 



2q ~ (37 - 2) = 3(2 - 7)E2 + (2 - 37)A^' + ,\,.2 "- 



Hence, it follows that for < 7 < 2/3, limr^oo 5^ = 1, and limr^oo V = Q. □ 

Theorem 2 (Global attractors for Ti{VIo)). For 2/3 < 7 < 4/3, any 
p € Ti{VIa) with n > 0, N'^ > N'^ and V < 1, has uj{p) = C. 

Proof. We note that in the subspace Ti{VIo) we have 

E+ = (, - 2)E+ - 2N^ - YTW^- (36) 
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Since q <2 (which can easily be checked) there exist a t\ such that < for 
all T > T\. Furthermore, in Tx{VIq) we also have 

'' = l+\,-l)V^ [(37-4) + 2E^]. (37) 

Hence, for 7 < 4/3, V will be monotonically decreasing for r > n, and 
linix^oo V = Q. This means that the solutions will approach the non-tilted 
subspace B{VIo) in the limit r 00. Thus after sufficiently long time, the 
monotone fmictions in Biyio) (given in [2]) will also be monotone along orbits 
in Tiyio). The point C is an isolated equilibrium point, thus we can use the 
monotone functions which implies uj{p) =C. □ 

There are strong reasons to believe that C is a global attractor in Ti{VI(i) 
for all < 7 < 2. As evidence for this is the existence of a monotonic function 
in Ti{VIq). Define a = E^. + S23, which in Ti{VIq) has the evolution equation 

a' = 2{q - 2)a - 4V3NNE23. (38) 

In Ti{VIo) we have the monotonically increasing function: 

7V2 + a ^, 4a(iV^ - jV2)(S+ + 1) ^ 
^' = WT^' {N^+a){N^+a) ^^^^ 

As T — > 00 there are thus three possibilities: 

<T, S+ + 1, or N'^-N'^ ^ 0. 

If cr ^ 0, the non-tilted analysis can again be applied. The case S_|_ = — 1 is 
unstable in the future so this cannot happen for general solutions. If A''^ — N"^ — > 
0, the tilted analysis of type II can be applied. All the late-time asymptotes of 
the type II case have an unstable direction into the interior of Ti{VIo) (which 
exactly corresponds to N'^ — N"^). Unfortunately, for the tilted type II model 
the local attractors have not been rigorously proven to be global attractors [8] . 

Apart from these two theorems we have not been able to show any global 
late time attractors for the various (sub)spaces. However, numerical analysis 
seems to imply that the local attractors also are global attractors. 



5 The initial singular regime 

Let us consider the initial singular regime; i.e. where r — > —00. This case is a 
lot more subtle than at late times due to the oscillatory behaviour of the system 
of equations as one approach r ^ —00. In fact, the tilted type II seems to have 
an initial oscillatory regime [8], and hence, one would expect a similar - if not 
more complex - behaviour for the tilted type VIq model. 

In the following we should emphasise that we consider a non-stiff fluid. For 
a stiff fluid (7 = 2) it has been pointed out that the Bianchi models allow 
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for a stable past attractor which would remove this chaotic behaviour into the 
past [28,29]. 

In the study of the initial singular regime it is useful to introduce the variable 
A, instead of N, as 

N = XN. (40) 

The equations for N and A are then 

N' = (g + 2S+ + 2\/3AS23)A^ (41) 
A' = 2\/3E23(l - A^). (42) 

In this case, the invariant subspaces A = ±1 correspond to T{II)'^. In this case, 
there are two Kasner circles, which differ by an orientation of frame. They are 
as follows.'^ 

1. /C±(//): Kasner "type 11" (0 < 7 < 2) 

+ Eig = 1, A = ±1, (7 = 2 
E_ = Ei3 = E12 = = O = 0. 

2. IC{VIo): Kasner "type VIq" (0 < 7 < 2) 

E2_ + E2_ = 1, A = 0, (7==2 
E23 = Ei3 = E12 - = = 0. 

Each of them also have extremely tilted Kasner sets and corresponding bi- 
furcations. These are discussed in detail in [30]. We will not discuss these here 
as they do not change our conclusion radically. However, bear in mind that 
these equilibrium points exist and that there are transitions between them. 

The behaviour near the initial depends on three types of heteroclinic orbits. 
These are 

1. Jr.: Frame rotations. 

2. 7jv±: Taub type II vacuum orbits. 

3. T\: Frame rotations between /C^(7/) and )C{VIo)- 

Let us consider, for illustration, the Kasner circle JC~^{II). The analysis for the 
Kasner circle K~{II) can be obtained from 1C~^{II) by using the map <j)2 (or 

h)- 

/C+(//) 

Here, the frame rotations, and Taub type II orbits are given by (when projected 

onto the (E_|-, E23)-plane): 

1. Three frame rotations: 
E+ = Ci 

E+ ± x/3E23 = C±. 
^With a slightly abuse of notation because they are really type I solutions. 
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Figure 2: The /C+(//) frame rotations projected onto the E23)-plane. Ar- 
rows are future-directed. 




Figure 3: Taub orbits, 7jv+, projected onto the (E+, E23)-plane. Arrows are 
future-directed. 
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Figure 4: Flow diagram for the transitions between the Kasner circles JC^{II) 
and IC{VIq) projected onto the (S23, A)-plane. Here, the constant r is given by 

r — y^l — E^. Arrows are future-directed. 

2. Taub type II vacuum orbits: 
E+ + 1 = C(S23 + V3). 

All of these orbits maps the Kasner circle onto itself: JC^{II) i~+ K.^{II). The 
frame rotations are illustrated in Fig|21 and the Taub orbits are illustrated in 
FigEI 

For the Kasner circle 1C{VIq) there are two frame rotations of a similar kind 
as for JC^{IT); namely the ones given by E+ ± \/3E_ = C±. These orbits also 
maps the Kasner circle onto itself: IC{VIo) i-^ IC{VIo). 

In addition to these homoclinic orbits there are homoclinic orbits between 
the three types of Kasner circles. They are given by the set of differential 
equations 

A' = 2%/3E23(l- A^) 
E'_ = -2\/3AE_E23 

E23 = 2\/3AE^ (43) 
1 = E^ + E^ + E23, E+ = constant. 

A flow diagram for this system projected onto the (A, E23)-plane is shown in 
Fig^ We can see that there are homoclinic orbits between JC^{II) and JC{VIq). 
However, these only form a set of measure zero. Most orbits connect JC^{II). 
Nonetheless, orbits can come arbitrary close to /C(V^/o) and hence, there might 
be frame rotations within the Kasner circle JCiVIo) before the orbit go back to 
/C±(J/). 
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From this fairly simple, but far from complete analysis, we believe that the 
initial behaviour of the tilted type VIq models is fairly complicated. There are 
infinite sequences of homoclinic orbits which consists of frame rotations and 
Taub vacuum orbits. These orbits are believed to be chaotic in general and 
thus we conjecture that the tilted type VIq has a chaotic behaviour to the past. 

6 Summary 

Here, for the first time, we have analysed the late-time behaviour of general 
tilted Bianchi type VIq universes. Our results are summarised in Table El We 
performed a local analysis and found all the future stable equilibrium points 
for various subclasses of tilted type VIq models as well as for the general tilted 
type VIq. In particular, we confirmed the observation in [1] on the existence 
of new self-similar solutions for 7 > 6/5. These solutions also proved to be 
important for the late-time behaviour; the Rosquist-Jantzen solutions are late- 
time attractors for a certain class of models with a two-component tilted fluid; 
and the line bifurcation at 7 = 6/5 is, in fact, the late-time attractor for general 
tilted type VIq solutions. For 7 > 6/5, the late-time attractor is an extremely 
tilted model.'* 

It is interesting to note that the general late-time attactors lie in the fixed- 
point-set of This means that in general the solutions are asymptotically 
(/<i-symmetric. More specifically, 0i-symmetric implies vi = 0, and V2 = '^l- 

At early times the analysis suggests chaotic behaviour as r —00. This 
behaviour seems to be a generic property of tilted Bianchi models [30] . So this 
behaviour was not very surprising, in particular considering the fact that the 
tilted type II model - which is part of the boundary of type VIq - was known 
to be chaotic [8]. 

Let us finish off with some comments about future research and how this 
work may relate to other Bianchi types. Firstly, the work [1] indicates that 
there may be some similarities between the more general type VI/i models and 
the type VIq. Hence, some of the features of the late-time behaviour found in 
this work may also appear in the type VI/j models. As goes for the class VHq 
model, which is given by N"^ > N'^, there is one obvious difference. The general 
non-tilted VHq model is not asymptotically self-similar [31,32]. This happens 
because the type VHq state space is not compact; there is no upper bound on 
N. In the terminology of [33] , the type Vllg model is extremely Weyl dominant 
while the type VIq model is Weyl-Ricci balanced at late times. Hence, we would 
expect a fairly different behaviour for the tilted type VIIq model at late times 
than that found here. Moreover, in this work the extremely tilted invariant sets 
have not been emphasised. For example, a two-fluid model where one fluid is 
extremely tilted is if particular interest. Further investigations in this case is 

*It should be noted that the 7 = 6/5 hne bifurcation looks remarkably similar to the 
Wainwright 7 = 10/9 line bifurcation of the exceptional model [27]. This makes us wonder if 
there are unknown line bifurcations along the entire line given by 7 = 2(3 + \/—h) / {5 + 3\/—h) 
for the more general type VI/i models [1]. 
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required. 



Nonetheless, there are quite a few unanswered questions regarding Bianchi 
models with a perfect fluid. This work has answered some of them. Hopefully, 
future work will answer more of them. 
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A Some simple proofs 
On bound dUD 

We start with assuming that the first two bounds in eq. (|27|l hold. Then, using 
the Schwarz inequality, we have 



Consider now the maximal value of the function F{X, Y) = AX^Y^ inside the 
unit disc 

+ Y^ < 1. 

One easily finds that the maximal value is at X^ — Y"^ = 1/2. Thus 

F{X,Y) < 1. 

Hence, by identifying X^ = and = E"! + 11^2 + ^13 ; obtain from eq. 
El 

4I]2_Ar2 + (Si2iV + SiaTV)^ + (SigiV + SiaiV)^ < 1. (45) 
Thus the last inequality in eq. (|27|l is redundant. 

Showing \S\ < 2S < 2 

We will consider the function 



where the matrix Y.ab is symmetric and trace- free, and c^Cq = 1. This implies 
that the maximal value of |5| occurs when c" is parallel to one of the eigenvectors 
of Eab. Thus, if Xi are the eigenvalues, we have 



< 4J:^_N^ + 2 (E?2 + ^13) (^^ + 




(44) 



\S\ = ISafcC^C^ 



SafccVl <max(|Ai|,|A2|,|A3|). 



(46) 
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The eigenvalues obey the relations 



S-^^ = Ai + A2 + A3 - 
Sai,!]"'' = A^ + A2 + A3 — 6S^. 

These two equations imply that the maximal value of an eigenvalue occur when 
Ai = ±2S, and A2 = A3 = =fS (or a permutation thereof). Hence, according to 
eqs.(@ni and we have 

\S\ < 2E < 2. (47) 

B Eigenvalues of Equilibrium points 

In this appendix we discuss some of the eigenvalues for the various equilibrium 
points. 

B.l Non-tilted 

1. T[I): FRW 

Ai,2,3,4,5 = -^(2 - 7), A6,7 = i(37 - 2). 

2. CS(1I): Colhns-Stewart type II (2/3 < 7 < 2) 

The essential unstable eigenvalue is for all equilibrium points in T{II): 

A7 = -4V3E23. 
Hence, since E23 < this point is unstable. 

3. C{VIo): Collins VIq (2/3 < 7 < 2) 

Ai,2 = -1(2 - 7) (1 ± , A3^4 = -1(2 - 7) (1 ± ' 

As = -§(2 -7), A6,7 = -|(6-57). 

Here, Ai,2,3,4 correspond to the non-tilted case, Xsfij are the eigenvalues 
in the vi, V2 and W3 directions, respectively. 

B.2 Intermediately tilted 

1. n{II): Hewitt's tilted type H (10/7 < 7 < 2) 
Unstable due to the eigenvalue 

A7 = -4\/3S]23. 
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2. £(//): Type II line bifurcation (7 = 14/9) 
Unstable due to the eigenvalue 

A7 = — 4a/3E23. 

3. 7^(F/o): Rosquist-Jantzen (6/5 < 7 < 3/2) 

Due to the complex character of these solutions, one has to part of the 
stability analysis numerically. Some eigenvalues are possible to find ana- 
lytically. The following seems to hold: 

Re(Ai,2,3,4) = -f(2-7), 
A5 + A6 = -i(2-7), Re(A5),Re(A6) < 0, 
A7 = i(57-6). 

Here, Ai_2,5,6 correspond to directions along the invariant subspace T2{VIo). 

4. C{VIo): Type VIo line bifurcation (7 = 6/5) 

Again we have to rely on some numerical analysis. However, analytic 
combined with numerics seem to indicate that 

Ai = 0, RC(A2.3,4,5) = -f , 

A6 + A7 = -f, Rc(A6),Rc(A7) <0. 
Here, Ai,2,3,4,5 correspond to directions along the invariant subspace F{VIo) 

.3 Extremely tilted 

1. £{II): Extremely tilted typo II (0 < 7 < 2) 
Unstable due to the eigenvalue 

A7 = -4\/3S23. 

2. £i{VIq): Extremely tilted type VIq (0 < 7 < 2) 

Ai_2 = -i (3±i\/T83) , A3,4 = -| (3 ± iVsV^S + SVS^ 

A5 = -!, A6 = -fc^, A7 = | (48) 
Here, Ai,2,5,6 correspond to directions along the invariant subspace T2{VIo)- 

3. £2{VIo): Extremely tilted type VIq (0 < 7 < 2) 

Ai = ^, A2,3 = -f(l±^^A9), 
A4,5 = -f (litVTT), A6,7 = -f (l±*\/l4 + 5\/2). (49) 
Here, Ai,2,3,4,5 correspond to directions along the invariant subspace F{VIo) 
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